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Abstract 

We experimentally demonstrate optical control of filamentation that occurs during propagation 
of intense, ultrashort laser pulses through crystals like barium fluoride and sapphire. Control is 
exercised by rotating the plane of polarization of the incident laser radiation and is demonstrated 
by directly visualizing filamentation in the bulk via six-photon absorption-induced fluorescence 
and, concomitantly, by probing the spectral and spatial properties of white light that is generated. 
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Filamentation is a spectacular nonlinear process that occurs when ultrashort pulses of 
intense light propagate through optical media. Its importance stems from the fact that 
optical energy is precisely localized within a filament. The ability to localize optical energy 
in bulk media opens tantalizing possibilities in areas like micromachining and laser surgery. 
Filamentation is also important in the macroscopic domain: in air it has been used to con- 
trol natural processes like lightening [1] and its study has been shown to be of utility in 
remote sensing |2j and in light detection and ranging pj • An important ramification of the 
light-matter interaction in condensed media is that the localization of optical energy may 
make possible localized modification of material properties like refractive index. This would 
open up diverse applications, like three-dimensional material modification and fabrication 
within bulk materials without surface damage; the feasibility of "optically fabricating" pho- 
tonic devices like resonators, beam-shapers, interferometers and amplifiers has already been 
demonstrated |4j. Periodic arrays of nanostructures have been reported in materials like 
fused silica p. Light-induced microexplosions within crystals have afforded applications in 
emerging areas like nanophotonics and plasmonics Such microexplosions arise from the 
plasma that is created when a small volume of the irradiated material (typically ~0.1 /im 3 ) 
absorbs energy from a femtosecond laser pulse of intensity ^10 14 W cm -2 . The resulting 
pressure (which can be in the TPa range) generates shock and refraction waves which, in 
turn, create nanovoids 

From a fundamental viewpoint, understanding the propagation of ultrashort light pulses 
in bulk media is complicated as both the temporal and spatial profiles of the incident pulse 
are altered due to linear and nonlinear effects like group velocity dispersion (GVD), lin- 
ear diffraction, self-phase modulation (SPM), self- focusing, multiphoton ionization (MPI), 
plasma defocusing and self steepening |7] . Filamentation and white light generation are two 
facets of propagation effects. White light generation is a manifestation of the dynamics due 
;o which additional frequency components are added to the incident spectrum. Recent work 
s| has shown that SPM is responsible for generating spectral broadening that is symmetric 
around the incident wavelength; at higher powers, other mechanisms like self-steepening and 
free-electron generation due to multiphoton excitation generate blue-side frequency compo- 
nents, resulting in asymmetric broadening. The physics governing filament formation in 
transparent media is yet to be fully understood. Currently, two broad modeling approaches 
are in vogue. One describes a filament as a soliton-like, self-guided beam wherein filamenta- 
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tion dynamics are governed by competing processes like the optical Kerr effect (that results 
in self- focusing of the incident light, thereby increasing its intensity) and multiphoton exci- 
tation (that generates electrons, creating a low-density plasma that, in turn, defocuses the 
beam) |9[. The alternative approach 10] visualizes filamentation in terms of optical energy 
being continuously absorbed and regenerated by subsequent focusing of different time slices 
of the incident laser pulse. 

Can some measure of control be exercised over filamentation dynamics? An answer in 
the affirmative would be of obvious importance, and constitutes the subject of this Letter. 
Introducing a periodic mesh creates periodic amplitude modulation of the transverse beam 
profile that leads to deterministic multiple filaments in liquids . Organization of regular 
filamentation patterns in air has been demonstrated by imposing strong fie ld g radients or 
phase distortions in the input beam profile ^^|, with a deformable mirror or using a 
double lens set-up Q]. Increasing the ellipticity of the input beam decreases the power 
that is required for multiple filamentation to occur For more than 30 years it has been 
believed that filamentation is initiated by random noise that breaks the axial symmetry of 
the incident beam [16]. Recent work has demonstrated that noise, by itself, does not play 
a significant role in the filamentation dynamics [17] and that deterministic vectorial effects, 
like polarization, may be of importance. We report here polarization-dependent propagation 
of ultrashort laser pulses in crystals like BaF 2 and sapphire. We have monitored the fila- 
mentation dynamics by imaging the fluorescence signal obtained upon six-photon absorption 
of incident laser light in the case of barium fluoride crystals. Such visualization yields di- 
rect evidence that polarization enables us to readily exercise spatial control of filamentation 
within the crystal. Polarization control is also demonstrated by concomitantly monitoring 
the spatial and spectral properties of white light that is generated. 

We used a Ti-sapphire laser yielding 36 fs pulses (820 nm wavelength) at 1 kHz repetition 
rate. We used incident powers of 10-300 /xJ, corresponding to 10-1000 times the critical power 
for self- focusing. A A/2 plate controlled the polarization angle of the light incident on each 
crystal. A 30 cm lens was used for loose focusing of the incident light. Crystal lengths were 
3 mm for sapphire and 15 mm for BaF 2 . The white light produced was imaged by a CCD 
camera and characterized by a fiber-coupled spectrometer (spanning the 400-1100 nm range). 
Filamentation was imaged using a digital camera and was analyzed using ImageJ software. 
Our crystals were randomly oriented and we rotated each around the light propagation axis 
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so as to maximize the generation of white light. This defined the axis around which our A/ 2 
plate was rotated. 

The BaF 2 crystal affords the possibility of directly visualizing the self-focusing/defocusing 
events in the filamentation process. This is due to six-photon absorption-induced emission 
(see the blue color filaments in Fig. 1). We recall that the band gap of BaF 2 is 9.1 eV 
and the incident photon energy is 1.5 eV. The images shown in Fig. 1 vividly demonstrate 
that the localization of laser energy within the crystal is amenable to ready spatial control: 
we spatially translate the filaments within the crystal by rotating the polarization vector, 
keeping the incident power constant. We quantify the polarization-controlled spatial motion 
in terms of z j, the distance from the entrance face of the crystal beyond which we observe 
the 6PA-induced emission; this is a robust indicator of the the filament-start position. z Q f 
has contributions from physical focusing by the external lens as well as from self-focusing 
due to n 2 , the intensity-dependent nonlinear refractive index of BaF 2 . We fixed the position 
of the external lens and rotated the A/2 plate. For data depicted in Fig. 1, the distance 
between the lens and crystal was 27 cm and the geometrical focus of the lens was located 
outside the crystal. With the A/2 plate at 0°, z /=3.5 mm (vertical white line in Fig. 1). As 
the incident polarization angle rotates there is gradual spatial translation of the filament- 
start position. At 45° the filament-start position lies well within the crystal, with z f ~5.9 
mm. The value of z Q f reverts to its original 3.5 mm when the A/ 2 plate is at 90°. The overall 
change in z Q f is significantly large, as much as 16% of the crystal length. 

Can the physical shift of z f be rationalized within the context of existing knowledge? To 
probe this, we take recourse to methods used in related numerical studies of light-induced 
damage in bulk media wherein z Q f has been related to experimental parameters [18] like 
focal length, size of the beam, the distance between the lens and the entrance face of the 
crystal, the incident laser power and P cr , the critical power for self focusing which, for 
isotropic materials, is (3.77A 2 )/(8 7rn n 2 ). We note that P cr depends inversely on n 2 which 
contributes to the net refractive index of a material: n = n G + n 2 |i?| 2 , where n is the linear 
refractive index. n 2 is further related to the material's third-order susceptibility tensor: n 2 
= (127r/n )x®. In case of cubic crystals like BaF 2 with space-group symmetry m3m, the 
effective value of is [3 : 



x (3 '(0 = z& + (A - 3A) |cos2( ! ) + 1 ' : 
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where 9 is the incident polarization angle. As 9 changes, the effective value of x^ a l so 
changes and this, in turn, alters n 2 . The n 2 value along the e-direction depends on the prop- 
agation direction while along the o-direction is independent of the direction of propagation. 
Thus, P cr is direction dependent. In order to account for the observed shifts in z c /, for fixed 
incident power, the value of P cr would be required to change sixfold! This is clearly an 
unrealistic proposition. We calculate that the change in n 2 for BaF 2 is only 0.43 xlO -16 cm 2 
W as 9 changes by 45°, resulting in only a 130% change in P cr . Other factors are clearly at 
play and need to be identified and explored. 

Nonlinear absorption effects will not affect the beam power at the onset of filamenta- 
tion because these will set in only after sufficient intensity has been attained. Indeed, the 
beam power may actually decrease with propagation due to GVD-induced pulse broadening. 
However, this effect does not significantly broaden the pulse as z Q f changes from 3.5 mm to 
5.9 mm. Other factors that we have considered include refractive index modification, color 
center generation, and induced birefringence. The loose focusing employed in our experi- 
ments rules out the possibility of modifying refractive index and generating color centers; 
melting [20] is also ruled out. However, we note that loose focusing might well give rise to 
cross-polarized wave generation wherein a single input wave is partly converted to a wave 
that is polarized perpendicular to the input plane polarization due to degenerate four wave 
mixing. Such a process has recently been shown to occur with high efficiency in BaF 2 j^lj 
and may lead to variation in the critical power for self-focusing, subsequently changing the 
onset of filamentation that is observed in Fig. 1. 

For uniaxial crystals like sapphire, the relation for effective x^ is considerably more 
complicated as there is now a direction-dependence for e-rays but not for o-rays |22|. The 
value of P cr will vary accordingly. In sapphire, where filamentation was not visualized due 
to the absence of six-photon-absorption-induced emission, we use white light spectra as the 
diagnostic. White light spectra from sapphire (see Fig. 2a) have two components: one due 
to SPM that arises from the Kerr nonlinearity (symmetric about the incident wavelength) 
and an asymmetric component that arises from MPI-induced free electron generation. The 
latter contribution is on the blue side of the spectrum and in all our measurements we 
saw a polarization-dependent change in this blue-side component. The Kerr nonlinearity 
that we observe depends, as noted above, on n 2 and x^ • The latter is a tensor with 81 
components. Depending on crystal type and symmetry the tensor may have many non-zero 
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components contributing to the effective which will result in change in the value of 
n 2 . Other experiments on isotropic materials have shown that at high intensities there are 
deviations from isotropy due to induced birefringence [23]: the induced birefringence makes 
any isotropic material behave like a uniaxial crystal. 

We also made simultaneous measurements of polarization-dependent white light spectra 
BaF 2 (Fig. 2b). The 820 nm light was blocked in this case; only the white light was 
recorded. Symmetric and asymmetric components are clearly visible. Recent studies have 
explored variations in the efficiency of white light production and of its spectral content 
as a function of incident laser energy [8] but, to our knowledge, studies of how the overall 
shape of the spectra change as a function of polarization angle have not been reported. Data 
depicted in Fig. 2c show that there is a major polarization-dependent change in the white 
light spectral distribution, specifically a reduction in the spectral width as the incident laser 
polarization nears 45°. As the blue part of the spectrum is predominantly due to plasma 
effects (free electron generation by MPI j^j]), the reduction in width sensitively maps the 
lessening of effective laser intensity within the medium which is driven by alterations in 
self-focusing conditions that are a consequence of changes in local refractive index. Over 
and above measurement of spectra, it is also instructive to directly visualize the white light 
pattern. Results are shown in Fig. 3 for BaF 2 irradiated at different polarization angles that 
demonstrate another manifestation of filament translation. When the filament-start is close 
to the crystal entrance the white light spot is large. As the half-wave plate is rotated, the 
filament-start position moves inside the crystal, reducing the white spot in the far field. We 
observed more than a two-fold reduction in the white light spot size when the polarization 
angle was rotated by 45°. At 0° the central white light zone is surrounded by concentric rings 
called conical emission; the rings depend on SPM in the radial direction. As the polarization 
angle is changed the ring contrast degrades and the spatial extent of the conical emission 
is reduced: in the radial direction only those frequencies that satisfy the phase matching 
condition form concentric rings. The diameter of the central white light spot is polarization 
dependent, showing a minimum at 45°. Similar results were observed for sapphire. We 
note that self-induced polarization in isotropic media has also been recently reported by 
Midorikawa and coworkers 

To place our polarization-dependent control results in perspective, we note that it has 
been possible to achieve a measure of control on filament morphology in air by forcing 
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aberrations in the incident light beam and by means of diaphragms l25j ] . Mysyrowicz 
and coworkers [|12j have shown that by varying the laser power the self-focusing distance in 
air can be altered. We also made measurements at various incident powers; Figure 4 shows 
how the onset of filament ation depends on laser power. At relatively low power we see a 
single filament whose start position lies well inside the crystal (4.6 mm from the front face). 
As the power increases, the filament-start position (z f) moves nearer the entrance face of 
the crystal. The most noteworthy feature is that in order to effect a change in filament-start 
position of ^3 mm, as much as a 5-fold increase in power is required. The superiority of 
our simple polarization method of controlling the position of the filament-start position is 
clearly demonstrated. Our findings have diverse ramifications in applied areas and in the 
basic sciences, particularly from the viewpoint of material modification in the bulk. 

Useful discussions with See-Leang Chin and Andre Bandrauk are acknowledged. The 
Homi Bhabha Fellowship Council is thanked for financially supporting one of us (JAD). 



[1] J. Kasparian et al, Science 301, 61 (2003). 
[2] J. C. Diels et al, Sci. Am. 277, 50 (1997). 

[3] F. Vidal et a/., IEEE Trans. Plasma Sci. 28, 418 (2000), and references therein. 

[4] C. B. Schaffer et al, Opt. Lett. 26, 93 (2001); A. M. Streltsov and N. F. Borrelli, Opt. Lett. 

26, 42 (2001); X. Wang et al, Appl. Opt. 43, 4571 (2004), R. Osellame et al, J. Opt. Soc. 

Am. B20, 1559 (2003); and references therein. 
[5] V. R. Bhardwaj et al, Phys. Rev. Lett. 96, 057404 (2006); C. Hnatovsky et a/., Opt. Lett. 

30, 1867 (2005); Y. Shimotsuma et al, Phys. Rev. Lett. 91, 247405 (2003). 
[6] S. Juodkazis et al, Phys. Rev. Lett. 96, 166101 (2006), Appl. Phys. Lett. 88, 201909 (2006). 
[7] R. R. Alfano and S. L. Shapiro, Phys. Rev. Lett. 24, 584 (1970), J. H. Marburger et al, Prog. 

Quant. Electron. 4, 35 (1975), N. Bloembergen, Opt. Commum. 8, 285 (1973), E. Yablonovitch 

and N. Bloembergen, Phys. Rev. Lett. 29, 907 (1972), A. Brodeur and S. L. Chin, J. Opt. 

Soc. Am. B 16, 637 (1999). 
[8] A. K. Dharmadhikari et al, Opt. Express 12, 695 (2004), Appl. Phys. B 80, 61 (2005), ibid. 

82, 575 (2006). 

[9] H. R. Lange et al, Opt. Lett. 23, 120 (1998), L. Berge and A. Couairon, Phys. Rev. Lett. 86, 



7 



1003 (2001). 

[10] A. Brodeur et al, Opt. Lett. 22, 304 (1997), M. Mlejnck et al, Opt. Lett. 23, 382 (1998). 
[11] H. Schroder and S. L. Chin, Opt. Commun. 234, 399 (2004), V. P. Kandidov et al, Appl. 

Phys. B. 80, 267 (2005). 
[12] G. Mechain et al, Phys. Rev. Lett. 93, 035003 (2004). 
[13] Z. Jin et al, Opt. Express 13, 10424 (2005). 
[14] G. Fibich et al, Opt. Express 14, 4946 (2006). 
[15] T. D. Grow and A. L. Gaeta, Opt. Express 13, 4594 (2005). 
[16] V. I. Bespalov and V. I. Talanov, JETP Lett. 3, 307 (1966). 

[17] G. Fibich and B. Dan, Physica D 157, 112 (2001), Opt. Lett. 26, 840 (2001), Phys. Rev. Lett. 
89, 013901 (2002). 

[18] M. R. Junnarkar, Opt. Commun. 195, 273 (2001), and references therein. 
[19] R. Adair et al, Phys. Rev. B 39, 3337 (1989). 
[20] V. Koubassov et al, Appl. Phys. A 79, 499 (2004). 

[21] A. Jullien et al, J. Opt. Soc. Am. B, 22, 2635 (2005), Opt. Lett. 30, 920 (2005). 

[22] Q. Guo and S. Chi, J. Opt. A 2, 5 (2000). 

[23] A. Owyoung et al, Phys. Rev. B 5, 628 (1972). 

[24] K. Midorikawa et al, Appl. Phys. Lett. 80, 923 (2002). 

[25] A. Dubietis et al, Opt. Lett. 29, 1126 (2004). 



8 



FIG. 1: Polarization-dependent control of the spatial location at which filamentation commences 
in a BaF2 crystal. The horizontal white line marked 15 mm indicates the spatial extent of the 
crystal. Vertical white lines are a guide to the shift in filament-start position with polarization 
angle. 

FIG. 2: Polarization-dependent white light spectra obtained upon irradiation of a) a 3 mm long 
sapphire crystal (50 /iJ) and b) a 15 mm long BaF2 crystal (140 /iJ). c) Polarization-dependent 
full width at half maximum (FWHM) of white light spectra from BaF2- 

FIG. 4: Variation of the onset of filamentation in BaF2 with incident laser power. The spatial 
control is expressed in terms of z G j, the distance from the entrance face of the crystal to the 
filament-start position. 



FIG. 3: Direct visualization of the size of the white light spot in BaF2 as a function of polarization 
angle (16 /iJ). Note also the angle-dependence exhibited by the conical emission. 
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